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The tunneling rate t v /h of a vortex between two pinning sites (of strength V separated by d) is 
computed using the Bogoliubov expansion of vortex wavefunctions overlap. For BCS vortices, tun- 
neling is suppressed beyond a few Fermi wavelengths. For Bose condensates, t v = V exp(— im s d 2 /2), 
where n 3 is the boson density. The analogy between vortex hopping in a superconducting film and 
2D electrons in a perpendicular magnetic field is exploited. We derive the variable range hopping 
temperature, below which vortex tunneling contributes to magneto-resistance. Using the 'Quantum 
Hall Insulator' analogy we argue that the Hall conductivity (rather than the inverse Hall resistivity) 
measures the effective carrier density in domains of mobile vortices. 

Details of vortex wavefunctions and overlap calculations, and a general derivation of the Magnus 
coefficient for any wavefunction on the sphere, are provided in appendices. 

PACS numbers: 03.75.Lm, 66.35.+a 



I. INTRODUCTION 

Mass and charge transport in superfluids, supercon- 
ductors and Bose Einstein Condensates (BEC) are gov- 
erned by mobility of vortices 0,0 ■ In two dimensions, 
vortex centers effectively interact as point charges in a 
perpendicular magnetic (Magnus) field. In the super- 
fluid phase, vortices are pinned at zero temperature by 
impurity potentials HQ. Their mobility, just below the 
pinning temperature scale, is dominated by thermally ac- 
tivated hopping 0, essentially following a classical Ar- 
rhenius law. 

At lower temperatures, quantum fluctuations may, in 
principle, admit tunneling of vortices under energy bar- 
riers, resulting in 'quantum flux creep' 0, El U ■ Ex- 
perimentally, vortex tunneling in superconducting films, 
manifests itself as low temperature magnetization relax- 
ation and non-activated, variable range hopping re- 
sistivity @. As magnetic field and disorder strength in- 
crease, vortex tunneling can turn into long range der- 
ealization. This amounts to a quantum phase transition 
from the superfluid into an insulating j 1 OL 1 1 lL 1 1 2| . or per- 
haps a Bose metal phase |l3j . 

A microscopic computation of vortex tunneling rates, 
has been an elusive theoretical goal. The semiclassical 
(instanton) approach EH, requires the determination of 
the 'vortex mass' in the presence of Magnus dynamics, 
and low energy superfluid phonons [l5L fl6l fl7l fl8l | . 

In the presence of short range, localized pinning po- 
tentials, it is simpler to compute (as we show below) the 
vortex tunnelin g ra te from the many-body wave- function 
overlap El El M E3 ■ 

In this paper we present detailed vortex overlap calcu- 
lations for the following systems: 

(i) The weakly interacting Bose Einstein Condensate 
(BEC). The tunneling rate between two localized pin- 
ning potentials of strength V, separated by distance d, is 




FIG. 1: Overlap of two antipodal vortex-antivortex pair states 
on a sphere, shifted by distance d, which are used in this pa- 
per to numerically evaluate the vortex tunneling rate. White 
circles depict directions of circulating currents. 



determined to be 

t v = V cxp (-C^n s d 2 + 0{l/n s e)) , (1) 

where n s , and £ are the boson number density and co- 
herence length respectively. Numerically C = 1 ± 0.02. 
The tunneling rate 0) applies to low density BEC's and 
their charged version, the 'Bosonic superconductor'. 

(ii) The BCS superconductor. Conventional super- 
conductors, with large electron density and core radius, 
have a tunneling rate which is suppressed by a factor of 
e -o.2k F d ^ wnere fo F i s th. e Fermi wavevector. Thus, one 
can conclude that BCS vortices cannot observably tunnel 
under barriers larger than a few Fermi wavelengths. 
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Our results point to the experimental regime where 
vortex tunneling in superconductors may have measur- 
able effects on magneto-transport. Basically, the super- 
fluid density should be much lower than the usual metal- 
lic electron density. This condition may be realized in 
cuprate (High T c ) films @, especially in the underdoped 
regim e l22l . and in highly disordered superconducting 
films |lllll2j. where phase fluctuations are important. 

The paper is organized as follows. 

Section[n]sets up the tunneling calculation by express- 
ing the two site vortex tunneling rate in terms of the 
pinning energy and ground states overlap. 

Section ITTT1 presents the overlap calculations of the in- 
teracting BEC (arriving at Eq. 0J), and the BCS super- 
conductor. 

Section IIVI presents the vortex transport theory and 
its relation to electrical conductivity in superconduc- 
tors. We use a quantum hopping Hamiltonian based on 
_ coupled to the QED field of superfluid phonons 
Following the theory of Ambegaokar, Halperin 
and Langer j2^| (AHL), we derive the variable range hop- 
ping exponent and temperature scale. The analogy to the 
'Quantum Hall Insulator' 0>HHIli3 is utilized to argue 
that the Hall conductivity of a bosonic superconductor is 
a robust measure of its boson number density n s . (This 
is in contrast to the commonly used assignment of the 
'Hall number' using the inverse Hall resistivity). 

We provide pedagogically instructive details in a series 
of appendices: the overlaps of vortex mean field state, in 
App. JXJ, the Bogoliubov theory, in App. |Bj), and Bo- 
goliubov de-Gennes equations, in App. (|B 3|) . and their 
matrix formulations in the spherical geometry, in App. 

In Appendix {Dj) we prove that the Magnus coefficient 
(adiabatic curvature) of an arbitrary wavefunction on the 
sphere is given by its average angular momentum density. 
This result, which is peripherally connected to the main 
subject of this paper, generalizes a previous proof con- 
necting the vortex Magnus action to the far field density 

m. 

We conclude in Section with a summary and a dis- 
cussion. 



II. VORTEX TUNNELING AND 
WAVEFUNCTION OVERLAP 

The relation between the vortex tunneling rate 
and wavefunction overlap follows the method of non- 
orthogonal ground states first introduced by Heitler and 
London pilj . 

Let us consider a homogeneous condensate described 
by a many-body interacting Hamiltonian H.q. Wc per- 
turb the system with two weak, and symmetrically sit- 
uated, pinning sites at positions x\ and at x 2 . The full 
Hamiltonian is 



where V, = An(ajj). The pinning potentials are repulsive 
for the bosons (A > 0) and are therefore attractive for 
the vortices. The separation between the pinning sites is 
d = \x-2 — Xi\. We assume for this exercise that a reflec- 
tion symmetry exists in Ti about a mirror plane between 
the pinning sites. The ground states of the two par- 
tial Hamiltonians satisfy (H + V 4 )|*°) = E \^°). Thus 
we can use the symmetric and antisymmetric superposi- 
tions as a variational ground states for the two symmetry 
sectors, 



l*° ± > 



|*?>±|*S) 



/2(i±|<*s|*g)|) 

Their corresponding energies are bounded by 
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(3) 



(4) 



The coherent tunnel splitting is defined as t v = (E+ — 
£P_)/2. To first order in the overlap, ( v I'i| , I / 2 ), wc obtain 
(see comment pflj): 



(* 1 |V 2 |*o)-(*o|*o)-(*?|V 2 |*°) 

7|(*S|*§>| , 



(5) 



where V = (*°|V 2 |*i) « Xn 1 (x 2 ). We neglect the term 
proportional to (\&] ) | V2 1 ^2) since it depends on the den- 
sity at the vortex center, which is assumed here to be 
small. 

Thus, Eq. J3J) establishes that the pinning potential 
supplies the 'attempt rate' of the tunneling. We shall 
see that the vortex wavef unctions' overlap decreases as a 
gaussian of their separation d. 



A. Overlap exponent and compressibility 

We briefly review the important result of Niu, Ao and 
Thouless (NAT) concerning vortex overlap. Con- 
sider a vortex centered at X, with core radius £. The 
Onsager-Feynman |3d | wavefunction is constructed from 
the uniform ground state tf'oj an d is therefore asymptot- 
ically correct for coordinates Xi far from the core center. 

$ x ~ n (exp (iftxi - X)) f(\ Xl - X\/0) % (6) 

i 

where f(y) — > 1 at y » 1. The overlap between two 
such wavef unctions, displaced by a distance d, is given 
by 

(*x\*X+d) = e -^exp(^ s d 2 f*fs{k)\ 

S(k) = — I d 2 x(%\Sn(0)Sn(x)\%) 
n s J 



H = H +V 1 +V 2 , 



(2) 



(7) 
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where W c is the overlap of the core area. </>{x) is the 
angle between x and the x axis. n s is the average density, 
Sn = n — n s , and k ~ 27r/£ is core wavevector cut-off. 

The structure factor S(k) is bounded by Bogoliubov's 
inequality [3l| 



S(k) < y/F(k) X (k). 



(8) 



For interacting bosons of mass m and velocity indepen- 
dent interactions, the equal-time correlator F is given by 



F(k) = ^{[n k ,[H,6n_ k }}} = ^-. (9) 

The compressibility \ is related to the sound velocity c s 
by 



lim x(k) — > ~. 

fc^o 2mc? 



(10) 



Inequality (JHJ) ensures that for the interacting BEC, the 
momentum integration in JJJ converges at low k, and 
therefore the vortex overlap exponent is finite in the ther- 
modynamic limit. 

We point out in Appendices JXJ and (|B 3(1 . that both 
the Gross-Pitaevskii (mean field) coherent state and the 
BCS vortex wavefunctions, suffer from a spurious overlap 
catastrophe, due to their unphysical, infinite compress- 
ibilities l32l. 



III. VORTEX OVERLAP CALCULATIONS 



We use boson coherent states IjAlll to set up a semi- 
classical expansion of the partition function, 

JdTjd 2 x(iip*d r Lp-n[^p* 7 ip}) 
o 



exp(-PE MF [$\) V[r,,r,*]e 



(13) 



The first exponential is the classical (mean field) energy, 
and the remaining path integral is over the fluctuation 
field r] = ip — 0. 

The classical field tp minimizes the variational energy 
(ip\H\ip). It solves the Gross-Pitaevskii (GP) equation 



(k{a) + v - /i + sl^)! 2 ) £0) 



= 



(14) 



With a weak pinning potential at the origin, and an ex- 
ternal magnetic field or rotation [3^ , a stable vortex so- 
lution can be found, whose approximate analytic form is 



<P 



uq r e 



(15) 



where (r, </>) are the polar coordinates of the vortex cen- 
ter. For numerical evaluation of the fluctuation spec- 
trum, the trial solution l|15f) must be improved upon by 
iterating Eq. {T3}. 



A. The Interacting BEC 

We consider a two-dimensional BEC with short range 
interactions, described by the second quantized Hamilto- 
nian 



n 



d 2 x 



(11) 

tjj^ (x) creates a boson of mass m and charge q at po- 
sition x\ (i is the chemical potential. The single parti- 
cle potential V(x) includes confining and vortex pinning 
contributions. The kinetic operator is 



ff(A) = J-(?V + *A 
lm \ % c 



(12) 



and A is a vector potential. In the case of bosons 
of charge q in a magnetic field B — B z and sub- 
ject to a rotation u) = w rot 2, the vector potential is 
A = {\B + ?nc<7 _1 cj ro t)i x x, and the single particle po- 
tential is shifted by AV(x) = (uj 1 - t/2c)(qB + mcLu ro t)x 2 . 

If the system is a uniform droplet of bulk density n s , 
the chemical potential gets pinned at u = gn s . The im- 
portant parameters of the condensate are the phonon 
velocity c s = y^gn s /m, and the coherence length is 
£ = h/mc s . 



B. BEC Vortex overlap 

In Appendix (0 the mean field vortex coherent state 
is shown to suffer from an overlap catastrophe. In 
Appendix JbJ, the Bogoliubov-fluctuations corrected 
ground state is given by 

|$) = A/'cxp(iV t <3V' t ) cxp(/7/. t ) |0) 
f(x) = <p(x)- [d 2 x'$*(x')Q(x',x), (16) 



where the operator Q is determined by solving Bogoli- 
ubov's equations. The density profile of \<P) is defined as 



5n(x) = ($|n(x) — raJ$), 



(17) 



and plotted in Fig. as a function of radial direction 
from the vortex core. 

The BEC vortex wavefunction overlap is given by 



|($i|<f> 2 )| = cxp (-Wo - Wi) 



(18) 



We note that the leading order Wo is proportional to the 
condensate density no, while the fluctuation correction 
W\ depends on d, but not on no- Wi(R) turns out 
to diverge logarithmically with system size R, which is 
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FIG. 2: The mean field condensate density profile no(r) cor- 
rected by the Bogoliubov fluctuations density Sn a (r) as a 
function of radial distance r from the vortex center, no is 
the far field asymptotic density. 




.5 

dy/n ii 



FIG. 3: The leading order overlap exponent for an antipodal 
vortex pair of a BEC on a sphere, displaced by a distance 
d. no is the mean field condensate density. Numerical data 
for Eq. l|B15fl are marked by crosses and circles representing 
different coherence lengths in the range £/R 6 [0.075.0.15]. 
The upper angular momenta cut-off was taken at i m ax = 60. 
The solid line represents the fit with C = 1 in Eg, 11911 . 



simply an artifact of the Gaussian approximation in | $) . 
At this sub-leading order, in l/(no£ 2 ), the compressibil- 
ity of this wavcfunction diverges with R. However, this 
need not concern us here, since we limit our calculation 
to the leading order which has a finite compressibility. 

In Appendix {BJ) we show that Wo of Eq. (|B16|) is 

dominated by the core region. A numerical calculation, 
described in Appendix Q, is used determine W(d, £). A 
fit to the numerical results for the antipodal vortex pair 
on a sphere of size R, depicted in Fig.QJ yields, 

W sphere (d, R) = Cirn d 2 (l + 0(£/R)) , C = 1 ± 0.02. 

(19) 

By halving the result for the vortex pair on a sphere 
(I19fl . an taking the large R limit, we obtain the leading 
order result for a single vortex in the plane 

W = \nn s d 2 + 0{i/R,5n), (20) 

which produces the result quoted in Eq.Q. 



C. BCS Vortex overlap 

In Appendix (|B 3|) the vortex wavefunction overlap for 
an s-wavc BCS superconductor is calculated by numeri- 
cally solving Bogoliubov-dc Gcnncs equations. Analysis 
of the far field contributions, combined with the numer- 



ical results yields 

|($i|$a)| = exp(-IF BCS ), 

H/BCS „ d 2^ Hm + k 2 pF{kFO y (21) 

where kp is the Fermi wavevector, £ is Pippard's coher- 
ence length and the vortex core size, and F is a dimcn- 
sionlcss function of magnitude 0.2. The overlap catastro- 
phe (logarithm of system size R) comes from the far field 
contributions, and is an artifact of the infinite compress- 
ibility of the BCS wavefunction. The important lesson 
from 1121 [I is derived from the core region contribution 
e -o.2k F d _ xhis suppression is not affected by the long 
wavelength phase fluctuations, but arises from the large 
density of electrons in the Fermi sphere which change 
their momenta when the vortex moves. We can safely 
conclude therefore, that in high superfluid density super- 
conductors, vortex tunneling cannot be observed beyond 
a distance of a few Fermi wavelengths. 

IV. VORTEX TRANSPORT THEORY 

Mobile vortices imply the destruction of the static su- 
perconducting order parameter. However when order 
parameter (phase) correlations are of sizeable range, it 
is still convenient to describe the dissipative and Hall 
transport in terms of dilute vortices, rather than in terms 
of a dense system of interacting bosons. In two dimen- 
sional BEC's, the neutral and charged BEC have analo- 
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gous transport equations, since shielding can be ignored. 
The number current, pressure and rotation frequency in 
the neutral BEC play the role of electrical current, volt- 
age and perpendicular magnetic field respectively, in the 
bosonic superconductor. For convenience, we shall dis- 
cuss the latter case, keeping a keen eye on possible ex- 
perimental ramifications. 



A. Vortex conductivity and electrical resistivity 

Steady dissipative vortex motion can be driven by a 
bias current. For bosons with charge q the current is j c = 
qn s v s , where v s is the supcrfluid velocity The vortex 
'charge' is the sign of its vorticity, Q v = ±1, and the 
vortex 'flux quantum' is unit of boson particle number 
= 1. The kinetic energy difference a vortex at site Ri 
and Rj is 

Eo(i) - E (j) = xz-iRi-Rj). (22) 

q 

The effective 'electric' (actually Magnus) field acting on 
that vortex is 



-Jc X Z. 



(23) 



The stable ground state of a charged superfluid (i.e. a 
superconductor) in a magnetic field B has a finite density 
of vortices given by 



n v = Bq/(hc), 



(24) 



where c is the speed of light. If the vortices have an 
average drift velocity V v the vortex current is given by 



. _BqQ 2L 
tic 



(25) 



The electromotive field (EMF) induced by the average 
vortex drift velocity is 

E = -c~ 1 V v x B = — — J v x z, (26) 

q 

where c is the speed of light. 

The 'vortex conductivity ' tensor is defined as o~ v , 



<Jv / j u v 
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(27) 



Using ( I23|) and l|25|) , the electric resistivity tensor p is 
directly related to the vortex conductivity by duality re- 
lations fill 



al\ (& y^x), 



- I - I o-« 

,q, 



(28) 



For example, consider unpinned vortices which can 
move in a Galilean invariant superfluid. The vortex Hall 
conductivity is given by analogy to a charged liquid of 
filling fraction j/, which is a xy = ve 2 /h, Setting e — > Q v , 
one obtains 



Qv ! 

h 



(29) 



which by using (|24|) and (|28fl yields (unsurprisingly), the 
classical electric Hall resistivity of a charged liquid with 
density n s : 



r class 



B 

n s qc 



(30) 



This result will be used later, in Section (|f V D|) . 



B. Vortex hopping hamiltonian 

Thus we arrive at non-interacting vortex hopping 
hamiltonian is described by an Anderson tight binding 
model in a strong magnetic field 



(31) 

where c] creates a vortex at a random pinning site 
with random energy a is the Magnus gauge field which 
satisfies 



V x a = hn^z 



(32) 



which gives rise to a Hall effect. The vortex hopping rate 
t v was precisely defined by JSJ. For a weakly interacting 
BEC it was shown to decay with inter-site separation d 
as a Gaussian 



t v (d) - V exp(- -n s d 2 ). 
The low energy vortex current operator is 



(33) 



J v (x) = T^^dijtvidij) (ct Ci - h.c.) S(x- Xij) , 

ij 

(34) 

where and Xij are the separation and midpoints of 
the pinning sites respectively. 

Vortex dynamics have been systematically derived for 
the two dimensional translationally invariant superfluid 
using an effective Quantum Electrodynamics (QED) the- 
ory |17tllSj- In the QED formulation, vortices are point 
'charges', moving in the presence of a 'transverse mag- 
netic field' (the condensate density), and interacting with 
'photons' (the Bogoliubov phonons). The photons give 
rise to a vortex self energy, which diverges logarithmi- 
cally at low frequencies. This diverging 'effective mass' 
however does not preclude quantum tunneling at finite 
timcscales I37TI . 
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Following the derivation of Ref. > we explicitly re- 
tain the low energy phonons by coupling them as a gauge 
field to the vortex current: 



H v ~ ph = J d 2 xJ v (x) ■ A ± (x) + he s \k\a f k a k , 



a _ \ " Ak x ( n s£ 

A ± - \ 9lfcl 



1/2 



z x k (a, 



2\k\i y-k^^-ki 



where creates a Bogoliubov phonon ('photon' in the 
QED language [l^) of wavevector k and frequency hc s k, 
and c s = h/(m£) is the speed of sound. The vortices are 
treated as hard-core particles, which like adsorbates on 
a surface, have a Fermi-Dirac occupation probability 



(>-">/ T + l) 



(36) 



We set the chemical potential /x to zero, and fix the av- 
erage density n v = X)i n */-^pin by the magnetic field as 
given by J23J. By inserting l(34|) in Eq. ( 135(1 for two sites 
separated by R, the two site conductance is given by the 
hopping theory ( pSSj l 

G"(d) = 7oe -™^ 2 e -2T(l e -l+l £ J l+l^-^D (37) 

where and 

2im v V 3 d l ^ phl 



7o 



/i 4 T 



■^(ci-ei), 



/•CO 

ft ph (w) = 3? / dte~ iwt (A x (x, t)A x (as, 0)} 

«/ DO 



4c s 



■(l + JV t (w/r)). 



(38) 



TZ ph is the bosons local dissipative response to the vortex 
motion, and Nj, is the Bose function. 



C. Vortex Variable Range Hopping 

The macroscopic vortex conductivity for the hopping 
hamiltonian (|31|) requires knowledge of the distribution 
of pinning site positions and energies V{{xi, £;}). Wc fo- 
cus our attention on individual vortex tunneling, in the 
regime of low vortex densities n v < < n p i n and small ran- 
dom fluctuations SV << V . Interaction effects between 
vortices are self consistently incorporated into e$. 

The density of states depends on both the pinning 
potential distribution, and the effects of vortex interac- 
tions. Here, we shall treat the low field, finite temper- 
ature regime where the mean hopping distance is much 
smaller than the vortex separation. The density of states 
is then reasonably modelled by 



Af(x, e) = S(x - Xi)S(e - e,)) = 



' ''pill 

W' 1 



(39) 



where n p i n is the pinning site density, and its energies are 
uniformly distributed in the interval (—5V/2,SV/2). 

The vortex conductivity maps onto Mott's variable 
range hopping (VRH) j3£| of charges in a random po- 
tential and a strong magnetic field. The gaussian decay 
of Eq. 1(33(1 is directly analogous to that of an electron 
in the lowest Landau level [40l Elj , where the 'Landau 
length' is A = \J2ixn s . 

At low enough temperatures T < SV, there are typ- 
ically many competing tunneling paths between pin- 



(35) ning sites separated by distances d >> A, 1/ 



by distances a >> A, i/^/n^T. As 
Shklovskii has shown |40t l4l| . in this regime multiple (vir- 
tual) tunneling processes play a crucial role. Their pri- 
mary effect is to replace the gaussian decay of the two-site 
tunneling rate by a linear decay, typical of an Anderson 
insulator. Here, wc must therefore replace t v (R) of 
by 



Ve 



(40) 



where £ is the linear localization length, and s is a numer- 
ical which depends on the details of P({a;,e}). Replac- 
ing ((33(1 by 1(40(1 . wc arrive at a two dimensional random 
resistor network, of the kind discussed by Ambegaokar, 
Halperin and Langer p3l | (AHL), with random conduc- 
tances given by 



7 exp 



2djj _ \ + \ej \ + \cj - ej 
I 2T 



(41) 



By AHL, the macroscopic conductance is given by the 
critical (lowest) conductance G% of the percolating subset 
of conductances which obey > G% . 

Taking the average number of bonds per site at per- 
colation to be v c (e.g. on the square lattice v c = 2), the 
percolating bonds all obey 



2d, 



i 2T 
which can be written as 

dij \ei\ + \ej 



< In( 7o /GS), (42) 



where 



2e E 



dmax = -\n(j a /G r v ), 



< 1, 



(43) 



Tln(7b/G$). 



(44) 



By l|39fl the density n con n of connected sites within 

| c z| _ c max 

is given by 



(45) 



SV 



The percolation condition on the number of connections 
per site is 



(46) 
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which implies the relation 

2 _ V C SV , 

e max a max — I 4 '/ 

"-pin 

Using (|44|l and l|47|l one obtains the value of the critical 
vortex conductance 

T = K5V(^) 2 , (48) 

where K = Av c j s 2 is a dimensionless factor of order unity. 
Using (|28|) we obtain the variable range hopping magne- 
toresistivity 

p**(B,T) = \ {n v {B))e-V°' T )K (49) 

p xx (B, T) exhibits vortex tunneling in two ways. First, 
the power of 1/3 in the exponential temperature depen- 
dence. Second, the VRH temperature scale To depends 
strongly on the ratio of vortex tunneling lcngthscales: the 
characteristic tunneling distance l/^/n p j n divided by the 
inter-boson separation \j ^fn~ a . 

At stronger fields (higher vortex density), long range 
vortex interactions are expected to modify the asymp- 
totic power of the hopping exponent 

D. Vortex Hall Resistivity 

1. Quantum Hall Insulator: review 

Early on, Holstein 0] nas studied the Hall effect of the 
hopping model (|31|l at low temperatures. He has shown 
the importance of three site tunneling interference, for 
producing a non zero Hall effect. 

Since then, several groups have extended that work to 
electrons in highly disordered two dimensional semicon- 
ductors in the presence of a strong magnetic field [24| ■ Al- 
though different approximation schemes were used, these 
groups have concluded that while p xx diverges at low 
frequency and temperature, liniT-»o,w-fO P xy ( UJ , T) < oo. 
Such behavior was dubbed 'Quantum Hall Insulator' (See 
comment |4*3|1. 

Experiments in Hall bars |44| , have found that the DC 
Hall resistivity has a much weaker temperature depen- 
dence than the resistivity on the insulator side of the 
field tuned metal-insulator transition. Refs. |25| have 
remarkably found the Hall resistance in the insulator to 
be quantized at h/e 2 v, at filling factor of v = 1/3. 

The apparent difference between the behavior of the 
Hall resistance versus the longitudinal resistance, can be 
explained by Kirchoff's transport theory of an inhomo- 
geneous resistor network |45^ , with widely vary ing ^ resis- 
tances. The Puddles Network Model (PNM) Jm, was 
introduced to explain the experiments of Refs. |25j. The 



PNM assumes a network of perfect Hall liquid puddles 
with conductivities a xy = ve 2 /h, a xx = 0, embedded 
in an insulating environment and connected by arbitrary 
large, classsical resistors. This model yields a quantized 
value of p xy = h/e 2 v : independent of p xx . 

While the PNM describes ohmic (incoherent) trans- 
port, quantum transport theory yields a different result. 
Using the Chalkcr-Coddington network to represents non 
interacting electrons in the lowest Landau level in the 
presence of smooth disorder, Ref. ^(| has numerically 
found that the Hall resistance at zero temperature ac- 
tually diverges with system size, similarly to the quan- 
tum induced localization of p xx . Therefore, a true Hall 
insulator phase for the Chalker-Coddington model has 
been ruled out. The conflicting results of classical and 
quantum transport theories, is related to the role of de- 
phasing. Inelastic scattering destroys localization, and 
prevents the divergence of p xy . 

2. Quantum Vortex Insulator 

The above discussion is directly relevant to the vortex 
hopping model 131fl . The vortices are essentially in an 
insulating state with possible domains of weaker super- 
conductivity where vortex mobility is higher. A diverging 
a xy at low temperatures may indicate long range coher- 
ent vortex transport, an interesting result in itself. 

Let us for now assume sufficient dephasing at the low 
temperature of experiments, due to vortex-phonon, or 
vortex fermion interactions. We can appeal to the Boltz- 
mann transport theory, and to the resistor network mod- 
els. This implies that the Hall conductance (not resis- 
tance) is determined by the Hall conductivity a xy of the 
most 'insulating' puddles. We do not know how to com- 
pute the distribution of a xy . However, by l|29(l . a xy mea- 
sures the effective carrier density n* in the most resistive 
domains, (i.e. the vortex liquid puddles): 

a xy {B) ~ n s qc/B (50) 

Furthermore, a detection of 'quantized' plateaux of 
o~ xy (B,n s ) may indicate locked-in charge density waves 
or topological ordering 01 i n the vortex-condensed do- 
mains. 

V. EXPERIMENTAL IMPLICATIONS AND 
DISCUSSION 

Cold atoms BEG. Vortices have been created in rotat- 
ing cold atomic gases 0]. One can imagine optically 
introducing localized pinning potentials and measuring 
the excitation spectrum. The lowest antisymmetric exci- 
tation could be compared to expression (Q, for different 
potential separations and boson densities. 

Cuprate superconductors. In thin cuprate films, time 
resolved magnetization relaxation Q , is a direct measure 
of the average vortex mobility. A variable range hopping 
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behavior of the magneto-resistance section is indicative 
of tunneling effects, as was shown in section <|1V C|) . For 
a 'bosonic superconductor' (coupled only to order pa- 
rameter phase fluctuations) , the characteristic resistivity 
given in Eq. (|49|l is 

<T T \npinj \T J \n pin J 

This expression has dubious applicability to high T c 
films. However, if one accepts a model of tighly bound 
hole pairs, with a low superfluid density, the vortices 
can primarily dissipate momentum to the low energy 
'nodal' fermions, and the core states near the vor- 
tex center [49l IHoj , Incorporating fermionic excitations 
amounts to adding a dissipative response to 7Z ph , given 
by 

U tol (uj) = (K^-) a fcrm (u>). (52) 

where er ferm (u;) is the fermions contribution to the AC 
conductivity. Crudely estimating the factors contribut- 
ing to the characteristic resistivity, we obtain 




The values of V and SV may be extracted from the re- 
sistance activation energy at higher temperatures. 

Hall Conductivity. Hall effect measurements in un- 
derdoped cuprates, have determined the Hall number as 
njj(T) = —B/(p xy ec) and found it to be of the same 
sign and magnitude as the hole do ping concentration 
away from the Mott insulator phase j5f]|. However, the 
'anomalous' strong temperature dependence of njj(T) 
has been used to distinguish the unconventional nature 
of the cuprates which differs from the much weaker tem- 
perature dependence of the Hall number in conventional 
metals. 

Our analysis, suggests that the Hall conductivity, 
Ea. H50p. rather than Hall resistivity should be used to de- 
fine a Hall number in the superconducting phase. In this 
regime, accessible by strong magnetic fields a xy is 
expected to be less temperature dependent and to charac- 
terize the Hall coefficient of metallic 'puddles' inside the 
superconductor, where vortices are locally delocalized by 
tunneling. 

Disordered Superconducting films. Highly disordered 
superconducting films \H\ are also likely candidates for 
observing vortex tunneling, since they effectively exhibit 
low superfluid density. We expect variable range hopping 
and a finite Hall conductivity near the supercondcutor- 
insulator transition, where vortices become delocalized. 
However, we refrain from quantitative estimates for these 
effects since a microscopic theory for strongly inhomoge- 
neous interacting fcrmion systems is beyond the scope of 
this paper. 



Periodic lattices. Optical lattices of cold bosons and 
Josephson junction arrays introduce the challenge of a 
strong periodic potential. The vortex hopping hamilto- 
nian (|3I|I can describe a periodic lattice of weak pinning 
potentials. One expects the lattice constant to play an 
important role in vortex mobility as it does in the Hof- 
stadter problem of a tight binding electron motion in a 
strong ma gnet ic field. Indeed, recent theoretical work 
has shown [47j that ground states degeneracies and vor- 
tex dynamics, depend on the boson filling per lattice site. 
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APPENDIX A: SPURIOUS OVERLAP 
CATASTROPHE IN MEAN FIELD THEORY 

This appendix provides pedagogical examples which 
reflect the limitations of mean field theory of superfluids 
of superconductors. In particular, we show that Bose 
coherent states and Bogoliubov de-Gennes wavefunctions 
have logarithmically divergent overlap exponents in the 
thermodynamic limit, as a consequence of their infinite 
compressibilities. 



1. Bose Coherent states 

Coherent states wavefunctions are often used as zeroth 
order approximations to ground states of Bosons, super- 
conductors, and quantum spin models with long range 
order. Here we show that when describing a vortex, co- 
herent states generically exhibit an overlap catastrophe, 
i.e. a logarithmically diverging exponent. This diver- 
gence is an artifact of the unphysical infinite compress- 
ibility exhibited by non- interacting bosons. 

Consider the coherent state with ip = y/noe 1 ^ a 
complex scalar field, defined by 

|^)=cxp Jd 2 x(tp(x)^(x) - tp*(x)i>(x)^)\0) . (Al) 

The overlap exponent between two translated vortex 
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FIG. 4: Overlap of two vortex coherent states as described 
by Eq. (S3). Black (red) arrows are vector representation of 
<pi (V2)- For asymptotic calculations of BEC and BCS cases, 
SCO Eas. HB16IB26J . the far field region is divided into uniform 
phase blocks labeled by I at distances r; and areas Ai. 



coherent states (|A1|) ~). centered at X = is then 



|<^2) I = exp(-W), with 



W 



cs 



d 2 x 



The far field integral (away from the core), where rii(x) ~ 
n . In this case, ip ~ ^fno exp(icj)), where (j> is the angle 
function relative to the vortex center (see Fig. g). 



no y d 2 a 

ro 



It 



\n 



cos 



-02) 



^nnod 2 ln(i?/r ), 



(A2) 



where ro is a near field cut-off. This calculation was 
first presented by Sonin |Tflj . in the context of vortex 
tunneling using a nonintcracting (product) condensate 
wavefunction. 

For the non interacting Bose condensate wavefunction, 
it is easy to verify that the structure factor is 



(6n(x)6n(x')) = S(x — x')no 



(A3) 



which yields finite zero momentum density fluctuations 
S cs (k) = 1. This is an artifact of the unphysical limi 
t of non interacting bosons, where the sound velocity 
vanishes and the compressibility is infinite. The wave- 
function overlap integral in Q diverges logarithmically 



with the lower momentum cut-off, signaling this orthog- 
onality catastrophe seen in l|A2[) . Recall, however, that 
non interacting bosons do not support stable vortices as 
their coherence length is infinite! 



APPENDIX B: BOGOLIUBOV FLUCTUATIONS 

The fluctuations 77 in (|13[) . are governed by the har- 
monic action 



where 



dr I dx (77* 77) 



(Bl) 



J = 



1 

-1 



Ho W 2N 
H (2) =1 I , (B2) 

,g<p* 2 H~ 



and Hq = K{A) + 2g \${x)\ 2 - /x 

The Hamiltonian is diagonalized by the canonical 
transformation, 



U V* 

V u* 



E 
E 



ith 



Here E 



&JS = J = SJS" 1 



(B3) 



(B4) 



m 'E n is the Bogoliubov spectrum [53|]. The 
spectrum and eigenopcrators are explicitly determined 
by solving the differential equations, 

H U n (x)+g^V n (x) = +E n U n (x) 



g<p* 2 U n (x) + H* V n (x) 



-E n V n (x) 7 (B5) 



1. Bogoliubov corrections in the far field 
approximation 

In a large area, asymptotically far away from the vortex 
core, we can solve the Bogoliubov equations i|B5|l using a 
constant order parameter, <p = -i/no e 1 ^ . is diagonal 
in Fourier space, and the matrix Q is given by 



QkM' 



= 2i0 



-- — (efe - M - E{e k ))8 k # 



' (B6) 



where e k = h 2 k 2 /2m and E(e) = y 1 e 2 + 2/xi. The Bo- 
goliubov fluctuations lead to an increase in the density. 
Defining 

Sn = ($\rfr)\$) = ±^ 2 (e fc ) 
k 

,2 



v\e) 



2V(e + A0 2 -M 2 (e + /x + y/e + f i) 2 - ]? 



(B7) 
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Changing variables to y = 1 + e//i, one obtains 



8n 



1 

7^ 

1 

4<2 



<'.v 



i 2y/y 2 -l (y+ vV-1 



(B8) 



The dimensionless parameter which controls the higher 
order terms in the saddle point expansion of Eq. (|13f) . 
is then Sn/no = l/(47mo£ 2 ), which serves as an effective 
'quantum disorder' coupling constant. 

For a uniform condensate, the leading order structure 
factor is given by 



■O(nor) 



2\-l 



(B9) 



where the order l/(no£ 2 ) corrections go beyond the lead- 
ing order approximation and hence ignored. An arti- 
fact of the wavefunction $, in i|lbj) is that it produces a 
non-vanishing contribution at zero momentum (i.e. an 
infinite compressibility) from the fluctuations correlator 
($|(?7) 4 |$). These are cancelled by self energy correc- 
tions due to cubic interactions gf*i]*(r]) 2 , which we will 
not calculate here. 

A quantum phase transition into another zero temper- 
ature phase (e.g. a solid) may be expected when Sn >n . 
Here we shall not explore the strong coupling regime. 



2. Full Vortex Bogoliubov Theory 

In Fig. [S] the numerical fluctuation spectrum about 
a vortex pair configuration on a sphere is plotted as a 
function of angular momentum. 

Eq. i|14|) has a zero mode [54| (E n = 0) corresponding 
to a global U(l) phase transformation cj) ^ (f> + S, given 
by 



U {x) = Atp(x) 

d 2 x(\u \ 2 -\v \ 2 ) 



Vo(x) = A* $*(x) 

. (BIO) 



This (unnormalizable) zero mode, which enforces charge 
conservation, is henceforth excluded from our numerical 
spectra. 

The Bogoliubov eigenoperators (quasiparticlcs) are 
given by 



d 2 x(u:(x)r,(x)-V:(x) V *(x; 



(Bll) 



To visualize the Bogoliubov fluctuations near a vortex, 
we plot in Figure 10 the tunneling density of states as 
defined by 

T(E,x) = 7rY / \Vn(x)\ 2 S(E-E n ), (B12) 



R 



m 



■ *■ »*4 ■ 
****** 



m 



FIG. 5: The exact Bogoliubov spectrum for a vortex pair con- 
figuration on the sphere, as a function of azimuthal quantum 
number m. The relevant parameters are £/R = 0.1, and a 
stabilizing rotation frequency of 0J Io t = 1.875 h/mR 2 . Note 
the non-degenerate zero mode at m — 1. Note the vortices 
induced distortion of the spectrum around m ~ in a broad 
energy range. 



which might prove interesting for inelastic scattering ex- 
periments of rotating condensates [55j . 

The Bogoliubov corrected ground state |$) is the a- 
vacuum, and thus satisfies a n 5?) = 0. This may be writ- 
ten in terms of the original bosons as 

3>) = A^cxp (\^Q^) exp(f^) |0) 
f(x) = <p(x)- fd 2 x' '<p*(x')Q(x' » (B13) 



Integration over coordinates is implied in (|BF 

1/2 

($|<I>) . The pair opera- 



and the 

normalization factor is N 
tor Q is 

Q{x, x') = (xKU^V^x') = Q{x\ x), (B14) 

in U and V is 



where elimination of the zero mode (jBlf 
implicitly assumed. 

We next compute the overlap of a vortex state with a 
displaced vortex. For two vortex states (|16f) centered at 
Xi, i=l,2. the operators Qi(x,x') are distinct and not 
translationally invariant. The magnitude of the wave- 
function overlap for Bogoliubov-corrected vortex states 
is 



($i|$ 2 )| 
W 



exp(-lT - Wi) 



W, 



(12) 



w, 



(11) 



Wi = - In 



D 



12 



IF, 



i -< 

-Qt i 



ft 



D, 



det" 1/2 (1 - Q*Qj) 



(B15) 
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T(E,r) 




FIG. 6: The Bogoliubov tunneling density of states of a BEC 
vortex (in arbitrary units) as a function of excitation energy 
E and radial distance from the vortex center r. £ is the co- 
herence length. 



We note that the leading order Wq is proportional to the 
condensate density no, while the fluctuation correction 
W\ depends on d,£,R, but not on no- Wi turns out to 
diverge logarithmically with system size. However, this 
need not concern us here, since we know that at the same 
order, the structure factor of $ is incomplete, as discussed 



after Eq. (|B9)l . 

We now show that in contrast to the coherent state 
overlap result (|A2|I . there is no logarithmic divergence, 
to leading order, in the exponent Wo(d,R). This can 
be shown analytically by separating contributions of the 
core and the far field regions, as depicted in Fig. 0] 

W° ~ f d 2 xw COIC (x)+J2Aiw taT ( Xl ) . (B16) 

core 

The far field integral is approximated by a sum over 
constant phase domains, at radii r/ 3> £, and of ar- 
eas Ai << rf, corresponding to the condensate field 
f = v^o ex P(^i,i)- We use the uniform solutions 1B6|) 
in each such domain, assuming that the block sizes are 
large enough: A x ^ £ 2 . Thus, 

p2ic/>i i 

Qu = -^-EQfce*^-*) 

' k 

= -e i2 ^(5(x-x')+0(Z/ri)) , (B17) 

which, by ifTflj). yields in each block I, and vortex config- 
urations i — 1,2, the constant functions 

fi,i(x) = (****•' -e-^ Jd 2 x'Q itl {x,x'^j 

= 2 y /r*e i +« +<D(t/R). (B18) 
Using (|B 1 5|) we obtain to zeroth order in £/r; the result 



Wfar^i) ^ no Re < 



a i2<t>2 



1 _ e i2(4>2-<t>{) \ e~' l(pl 



> = 4n Re 



„i(02-0l) 

e l 

1 _l_ e i(4>2-4>i) 2 



= , 
(B19) 



where we have suppressed the block index I. Unlike in 
the coherent states overlap exponent which exhibits a 
logarithmic divergence (|A2() , Wq is perfectly finite in the 
large system limit 

CO 

Aiw\ ar = fd\(d- V(f>) 2 x O (£/r) < oo (B20) 
1 ro 

which is in agreement with NAT (7J and (|B9|I . 

3. Vortex Overlap in BCS States 

It is perhaps little appreciated in the literature that 
vortex wavefunctions of BCS theory of superconductors, 



as derived by leading order Bogoliubov de-Gennes equa- 
tions, suffer from the same overlap orthogonality catas- 
trophe and infinite compressibility as the Bose coherent 
states described above. 

We shall demonstrate this point by considering a 
'generic' BCS superconductor, which is described by 
a microscopic Hamiltonian of attractively interacting 
fermions: 

n =Jd 2 x {^] (k + v) ^ ~ f v44Vv^ s } 

where ijj s (x) creates an electron of mass m, charge q and 
spin s =T,1 at position x, and the operators of kinetic 
(K) and potential (V) energy are the same as defined for 
bosons in l|12lll|l . Summation over repeated spin indices 
s, s' is assumed. The complex superconducting order pa- 
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rameter is f{x) = (ip^(x) ijj^(x)). At long distances from 
the edges or vortex cores, tp minimizes the Ginzburg- 
Landau energy, i.e. it satisfies Gross Pitaevskii equation 
(|14fl with pair mass 2m and charge q = 2e. Its magni- 
tude is given by the BCS gap parameter, i.e. A = g\ip\. 
As for the BEC, in the presence of a weak magnetic field, 
a quantized vortex solution minimizes the mean field en- 
ergy, and its core size is given by the coherence length 
£ ~ Hvp/ttAq, where vf is the Fermi velocity Q. 

The Bogoliubov de-Gennes (BdG) equations for the 
superconductor are 

H U n (x) +g<pV n (x) = E n U n {x) 
gip*V n (x) - H*U n (x) = E n V n (x), 

where Hq = K — fj, and /i = h 2 k F /2m, along with the self 
consistency condition, 



(B21) 



The self consistency determines the detailed profile of 
|<^(t)| in the core region. 



We first calculate the far field contribution, and see 
that it exhibits from a similar logarithmic divergence as 
the boson mean field wavefunctions. Asymptotically far 
from the cores, whose sizes are of given by £, one can 
diagonalize l|B21f) using a constant order parameter (fi = 
Ae^. The solution of jjESU) yields 



vi = 



A 2 



tt + A 2 



(B24) 



(B25) 



with (fc = h 2 (k 2 — kp)/2m. Factorizing the determinant 
into blocks, as shown in Fig. 21 wc obtain 



W 



,,far 



W core + ^ Ai W'i 
ri>r 



.far 



ulvd 



w 



[(l-COS (4>l.l - 02,i)] JJ^2 U k v k 

±tt n cS d 2 \n(R/r Q ) + W cole , (B26) 



d 2 k 




E, 




m 



FIG. 7: The Bogoliubov-de Gennes spectrum for a vortex 
pair configuration on the sphere, as a function of azimuthal 
quantum number m. The relevant parameters are the BCS 
gap A, and the Fermi energy £f- Note the branch of (doubly 
degenerate) low energy Caroli-de Gennes-Matricon core states 
at positive angular momenta. 



where n e ff = h-p/Awt; < n e . 

The divergence of the first term with log(i?) arises from 
the far-field contributions. Again, it is an artifact of the 
BCS wavefunction, which is a mean field state where the 
condensate phase does not fluctuate. As for Bose coher- 
ent states, the BCS wavefunction also has a non vanishing 
structure factor at zero momentum, and thus an infinite 
compressibility at zero temperature. This divergence is 
cancelled by including (RPA) phase fluctuations, which 
restores a finite zero temperature compressibility given 
by the density of states at the Fermi energy. 

Here we are interested in the core contribution Wc 0re 
of Eq. (|B26p . This requires a full diagonalization of 
the BdG equations. For the antipodal vortex pair on 
the sphere, the details of the computation are found in 
Appendix JCJ). A fit to the numerically obtained values 
of W in Eq. 1|B23(1 yields the asymptotic expressions at 
large R/d: 



W 



BCS 



k 2 F F(k F £) ) +0(d 4 ), (B27) 



Using the solutions of (|B21|) the vortex ground state is 
given by 



|$) = A/"exp(V[Q^) |0) 
Q = (xKU^V^x 1 ) , 



(B22) 



where Af is the normalization. The normalized overlap 
of two BCS vortex states displaced by d is given by 



-W(d) 



det 



(u\u 2 + V?V 2 



(B23) 



where F is a dimensionless function of order 0.2, of the 
scaling variable as demonstrated by the collapse of 

the numerical data for W in Fig. ||HJ). 

The first i?-dependent term is the diverging far field 
contributions discussed earlier, which is an artifact of the 
infinite compressibility of the BCS wavefunction. 

The second term is the core contribution which is in 
fact very large: it goes as (kfd) 2 . Thus we find that 
the core contribution effectively suppresses energy con- 
serving tunneling between sites separated by more than 
a few Fermi wavelengths. In the regime of large coher- 
ence length relative to the mean free path, vortices could 
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FIG. 8: The function F(k F ,£) denned in Ea. l|B27jl . which 
yields the coefficient of the core contribution to the BCS vor- 
tex overlap on the sphere. k F is the Fermi wavevector and 
£ is Pippard's coherence length. The numerical data in the 
specified ranges of k F R and (,/R, collapse to a function of the 
product ^k F . 



be dissipatively dragged by the bias current by exciting 
low energy core states [f5(| ) . This is essence the source of 
friction in Bardeen-Stephen flux creep theory |Ht| . This 
'factional' motion involves impurity scattering inside the 
vortex cores [H^, 59], which we do not treat in this paper. 



(1B15|) into a computationally convenient form 



x (JJ\ (<p 2 - Vl ) - V] (<p* 2 - ri)) 



W 1 

s 



|lndet[5ia] 
u\ U 2 - Vi v 2 



(C2) 

For the vortex pair field tp given by i|Cl|> the BdG equa- 
tion possesses axial symmetry which allows m = fh n to 
be a good quantum number, and the BdG eigenvectors 
are 



u l,m 
y l,m 



JJ n S 

^ I u m,m 7l i 

EW- 



for these coefficients is 

jjn \ / H N 

-V n ) ~ \-A* 
H N {m) - h ' 

VlV,m 

A w (rn) 



A 

-H N 



jjn 

-V r ' 



(C3) 



(C4) 



nR ?l(l + 1)5,,, + V w , m 

(l,m\V pin (9) + 2g\^\ 2 - fi,\l',m) 

g(l,m\{p 2 \l',m-2) (C5) 

where ()C8ll can be used for the precise numerical evalu- 
ation of the matrix elements. 



APPENDIX C: CALCULATIONS ON THE 
SPHERE 

An antipodal vortex pair field can be expanded as 

oo 

p(M) = X>>li(M) (ci) 
i=i 

where Yi m (9, </>) are normalized spherical harmonics. Eq. 
(|C1|) is an eigenfunction of L z = —ihd^,. The coherence 
length £ determines the core sizes at the north and south 
poles, and thus the the decay rate of the tpi with I. 

The angular momentum representation, with a cut- 
off at Z max 3> (-??/£), reduces the Bogoliubov equations 
(|B5IB21ll to finite matrix diagonalizations. We can also 
translate the vortex wavefunctions using SU(2) rotation 
D Lm> matrices. 



2. Boson Ground State Overlap 

Having determined C/,",VJ", and recognizing that the 
transformation between (fo,Uo.Vo and <pi,Ui.Vi simply 
involves an 0(3) rotation of the z axis by an angle 9, the 
overlap (|C2|I is given by 

det- 1/2 [S1 exp [^IS^J^] 



($i|$ 2 ) 

ll,n 
Jn,l 



D ln n -2,-l) 



= dnYMe,4)<p(9,(/>) 



(C6) 



ft is easy to verify that for the limit UJ 1 = 5 n j, V, n = 0, 
J6|) reduces to the result for free particles. 



1. BEC with Antipodal Vortex Pair 

For the BEC vortex pair, some straightforward but 
tedious algebra, can bring the overlap exponent in Eq. 



3. Fermions on a Sphere 

A polar vortex pair field described by Eq. (|C1|) , defines 
the pairing order parameter 



4 = £(M) 



(C7) 
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by normalizing it such that -v/no = Aq, In the spherical 
harmonic basis |Z, m) the Hamiltonian (|B21|) simplifies 
greatly. The Laplacian is proportional to the diagonal 
operator L 2 , and non diagonal matrix elements of func- 
tions F(0, </>) can be computed using 3 7 Racah coefficients 
163: 



(l,m\F \l',m- M) 

OO 

L,M 

x(-iy 



(2Z + 1)(2£ + 1)(2Z' + 1) 

47T 



1/2 



I L V \( I L V 
-m M m- M \ 



F LM = fdnY£ M (fl) F(O) 



(C8) 



Due to axial symmetry, to is a good quantum number, 
which is to say that a function m n is defined such that 



v; 



^ l v m,m rl ; 
v l> u m' ,m„- 



(C9) 



The matrix BdG equation for these coefficients is 
V" ) ~ \ At -H N ) \ V r 

An' (to) = (Z,m|£|Z',m- l) (CIO) 



The overlap of two vortex pair states relatively rotated 
by 9 is 



($i|$ 2 ) 



det 

nil' 



E { u ? Dl m 



(Cll) 



+V, n D 



I m n — l,m n f — 11 



v,' 



where D l mm , (9) is the orthogonal rotation matrix. 



APPENDIX D: THE MAGNUS ACTION 

The calculation of the Berry phase for the motion of 
a vortex wavefunction can be done using the method of 
Arovas, Schrieffer, and Wilczek, [6l|. originally for quan- 
tum Hall effect quasiparticles, and later applied mutatis 
mutandis to supcrfluid vortices by Haldane and Wu • 
Here we show how the Berry phase is calculated simply 
and exactly on a spherical geometry. 

The spherical geometry yields many advantages for 
vortex wavefunctions. The geometric phase of a moving 
vortex is tricky to evaluate for superfluids and supercon- 
ductors on a finite plane, since it sensitively depends on 
the boundary conditions. On the sphere, there are no 
boundaries to worry about. The translations are imple- 
mented by 0(3) rotations, whose generators do not com- 
mute. Hence we can show that the Magnus density of a 



general many-body wavefunction is simply given by the 
expectation value of angular momentum density. When 
applied to an antipodal vortex pair state, this results 
agrees with Thouless, Ao and Niu's conclusion j23| that 
the Magnus density of a single vortex is just the super- 
fluid density away from the vortex cores. 

Consider a general many body wavefunction l^o) de- 
fined on a sphere of radius R, and calculate the Berry 
phase acquired by an infinitesimal loop of area A = 
R59 ■ R8(f) in parameter space, as depicted by Fig. |UJ 




FIG. 9: The orbit of a vortex pair on the sphere which is used 
to calculate the vortex Berry phase and Magnus action. Mo- 
tion of the antipodal vortex pair is achieved by four successive 
rotations of the north pole, covering a solid angle of x <j>. 



The infinitesimal loop is divided into four segments 
(i, i + 1) such that the loop Berry phase is 



(Dl) 



For a general wavefunction, the loop can be defined by 
a succession of small 0(3) rotations of a vector passing 
through a special point on the sphere. For our purpose, 
we pick that point to be the vortex center. For a gen- 
eral wavefunction, however, the choice of special point is 
somewhat arbitrary. Inversely, for a given sequence of 
small rotations, each point on the sphere executes a an 
orbit. The special points on the sphere which character- 
ize the sequence are the two antipodal points whose loops 
have a maximal area, (see Fig. |5J). Without loss of gen- 
erality, we choose to place these points at the north and 
south poles, and and apply the corresponding sequence 
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of rotations: 

|Vh) = e^/^o) 

|V> 2 ) = e iL * s<l>/H e iL * 50/h \rj} Q ) 

Uj)^ = e -iL y Se/h e iL x 5<j>/h e iL y 89/h 

|^4) = |Vo) (D2) 

The overlaps uto quadratic order in 50 and 5(j> is found 
to be 

e iSu/h w l-i666<j>(L z }/h + ... (D3) 

where the expectation values are in the state \ipo)- 
Generalizing to a complete path, we obtain 

S M = - [dt(Mt)\L\Mt))~, (D4) 

where the vector <1> is the rate of change of solid angle, 
instantaneously directed toward the vortex core. 

This Magnus action applies to any wavefunction, re- 
gardless of its correlations, (e.g. superfluid, Fermi liquid 



or solid). For an antipodal vortex pair in a superfluid, 
the angular momentum density is given by h times the 
condensate number density no, if one assumes that the 
thermal excitations (normal component) carry no angu- 
lar momentum, representing the vortex center coordi- 
nates by (X,Y), yields the Magnus Lagrangian 

C M = 2Trhn Q XY (D5) 

which resembles the effect of a uniform magnetic field in 
the z direction. 

The fluctuations correction of the Bogoliubov cor- 
rected wavefunction (|16f) is given by 

Cm = 2irh (n + 5n) XY (D6) 

Interestingly, although Sl z (x) ^ Sn(x), the system aver- 
aged quantities are found to be numerically equal Sn = 
81,. 
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